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Abstract 

We reanalyse the production of pp pairs in proton-nucleus and nucleus-nucleus 
collisions employing novel elementary cross sections for baryon-baryon and pion- 
baryon production channels based on a boson-exchange model. In contrast to 
previous transport studies performed in the literature the secondary pion induced 
channel is found to be most important in both p+A and A+A collisions at sub- 
threshold energies. A detailed comparison with the experimental data available 
indicates that sizeable attractive p potentials in the order of - 100 to - 150 MeV 
at normal nuclear matter density are needed to reproduce the size and shape of 
the experimental spectra. 
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1 Introduction 



The production of particles at energies below the free nucleon-nucleon threshold ('sub- 
threshold production') is one of the most promising observable for hadron selfenergies 
at high nuclear densities since the particles are produced predominantly during the 
compressed stage 0, ^ (3J. Antiproton production at energies of a few AGeV is 
the most extreme subthreshold production process and has been observed in proton- 
nucleus collisions already a few decades ago [Q, ||, [/J. Experiments at the JINR || and at 
the BEVALAC ||, |10| have provided, furthermore, first measurements of subthreshold 
antiproton production in nucleus-nucleus collisions followed by measurements at KEK 



nj] and GSI [L3|] with new detector setups. Various descriptions for these data have 



been proposed. Based on thermal models it has been suggested that the antiproton 
yield contains large contributions from AN — > p + X, A A — > p + X and pp — > pN 



production mechanisms |L4|, y| ]lq] . Other models have attempted to explain these 



data in terms of multiparticle interactions [17]. 

Nowadays, the relative strength of the various production channels as well as possi- 
ble in-medium effects are most effectively controlled by means of transport approaches 



IS, 19, EQ, EH E2L E3f. First results of a fully relativistic transport calculation for an- 



tiproton production including p annihilation as well as the change of the quasi-particle 



properties in the medium have been reported in [21 1. There it was found that according 



to the reduced nucleon mass in the medium the threshold for p-production is shifted to 
lower energy and the antiproton cross section prior to annihilation becomes enhanced 
substantially as compared to a relativistic cascade calculation where no in-medium ef- 
fects are incorporated. A variety of transport studies have been performed since then - 
most of them using the parametrization of the elementary production cross section as 



proposed by Batko et al. |19[ - and confirmed the necessity of attractive p potentials, 



however, with some debate about its actual magnitude |L8|, [20], |23| . 

More recently, the question about the elementary p production amplitude has been 
addressed by Lykasov et al. in an effective boson exchange model |2~4fl. Here the 



production cross section close to threshold follows accurately the 4-body phase-space 
constraints, but is much lower than the commonly adopted parametrization from [|l9fl . 
Such a reduced production cross section in baryon-baryon collisions would demand for 
even stronger attractive p potentials if further reactions channels might be neglected. 
However, as has been found in studies of K + and K~ production |25j |26[ the sec- 
ondary pion induced channel might dominate very well such that also the nN — > Npp 
production channel has to be considered. 

The paper is outlined as follows: In Section 2 we briefly recall the main ingredients 
of the relativistic transport approach used for p production before [p3[] . Section 3 is 



devoted to the calculation of the reaction nN — > Npp in a boson exchange model and a 
presentation of the elementary exclusive and inclusive cross sections. Section 4 contains 
a detailed comparison of the calculated spectra with experimental data for p+A and 
A+A collisions with the aim to obtain some closer constraints on the in-medium p 
potential, while a summary on the p selfenergies concludes the paper in Section 5. 
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2 Transport approach for antiproton production and 
propagation 

The phase-space distribution function for the antiprotons fp(x,pp) is assumed to follow 
an equation of motion equivalent to the transport equation for baryons [^7| . 



1^ 



^ + (n»+m;(a>;))9{; f 



fp(x, rip) 



with 

and the mass-shell constraint 
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f p (x,U p ) = 0. 



(2) 



(3) 



Since the momentum dependence of the scalar and vector selfenergy of the antiproton 
are unknown, we assume these quantities to be independent on momentum in Eq. (1), 
however, to be of different strength as compared to the nucleons, 



Uf(x) 
Ul{x) 



-9 P S o H {x) 
9v vh(x), 



(4) 



where a and u are the scalar and vector meson fields in the nuclear Lagrangian (cf. 
Ref. |p3| ). Thus the effective mass and the effective momentum of antiprotons are 
given by 



m*= = m + U$(x) 



IP 



rt-U^x) =f lx -glu«{x). 



m ~ 9%cr H {x) 



(5) 
(6) 



According to the arguments given before the coupling constants g v s and g% are treated 
as free parameters that should be fixed in comparison to data. 



2.1 The collision integral 

The collision term I^ oll (r.h.s. of eq. (|I|)) includes i) a term /^ ast , describing the elastic 
antiproton-baryon scattering, ii) a term Ip rod , describing the antiproton production 
by baryon-baryon and meson-baryon collisions, and iii) a term I^ bs responsible for 
in-medium p-absorption on baryons. The production of pp pairs by meson-meson 
collisions can be neglected in the energy regime of interest here since the invariant 
energies available are too small. l v e i ast describes elastic baryon-antiproton scattering 
as well as elastic antiproton-antiproton scattering. While this part of the collision 
integral l v coll can be formulated similarly to the collision term describing bayon-baryon 
scattering the other terms represent extensions to the conventional collision integral. 

The basis for the description of the p production by baryon-baryon collisions is the 
process 

5 + 5^ + ^ + ^ + ^ = 1 + 2-^ + 3 + 4 + 5, (7) 
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for which the corresponding covariant collision integral reads 



Iprod( X i 



m? mA rrin m A mi m; 



d 3 nid 3 n2d 3 n 3 d 3 n 4 d 3 n5- 1 ° 1 : ' " 



(2tt) ii 7 1 * d 4 n?n°n°n°n°n° 
w(n?, n£ | njf , n£, n£)5 4 (rf + ^ - ^ - ^ - ^ - njf) 
{/(x, no/(x, n 2 )(i - /(x, n 3 ))(i - f{x, n 4 ))(i - /(x, n 5 ))} , (8) 

where IV is the transition probability for the reaction Hi + Il 2 — > n 3 + n 4 + n 5 + lip 
in terms of momentum coordinates. We have omitted the Pauli-blocking factor for 
the antiproton in the final state because the number of antiprotons created during a 
heavy-ion collision in the subthreshold energy regime is negligible. For the same reason 
we neglect the effects of the reaction (0) on the phase-space distribution function of 
the baryons. The production term from irN collisions has a similar structure as Eq. (|8|) 
except for the fact that particle 2 denotes a pion and the final state is of 3-body type 
(with 2 nucleons only) since the 4-momentum of the impinging pion is used entirely 
for the production of the pp pair. 

In the collision integral I^ bs we do not treat all possible annihilation reactions sep- 
arately but sum up all channels in the inclusive annihilation reaction 

B+p^X, (9) 

where X denotes all possible final states (essentially pions) of the baryon-antiproton 
annihilation. The corresponding energy and momentum conservation reads 

p"b+p" p =Px, (10) 

where and pp denote the 4-momenta of the baryon and the antiproton, respectively, 
and Px denotes the sum of the 4-momenta of all particles in the final state of the 
annihilation reaction. Reaction @ then leads to the collision term 



(2tt) 3 J 1 A n? n°- 
w{ut, n£ | n x )5 A w +p$- Px )f(x, no/Or, n p -) (11) 

with W(ni,n^ I 11^) denoting the transition probability for the reaction (|9|). Inte- 
grating flnp over d 4 Hx we obtain 

Here the integration over d 4 Ux implies in addition to the integration over all final 
momentum states of a particular reaction a summation over all possible annihilation 
channels. Due to this fact , 11^) denotes the probability of an antiproton with 

effective momentum to annihilate with a baryon with effective momentum 
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2.2 Numerical implementation 



Since the production probability for antiprotons is very small the average time evolu- 
tion of the nucleus-nucleus is not affected and it is justified to treat the p-production 
perturbatively p], [|. Within this approach the p invariant differential multiplicity 
is obtained by summing incoherently over all baryon-baryon and meson-baryon colli- 
sions and integrating over all residual degrees of freedom. Assuming the antiproton 
production to take place via reactions of the type 

5 + 5^p + p + iV + iV=l + 2^p+3 + 4 + 5 (13) 



7r + 5^p + p + iV = l + 2^p+3 + 4 (14) 

(B stands for either nucleon or A) the invariant multiplicity as a function of the impact 
parameter can be written as 



BBcoll 



;i-/(x,n^))(i-/(x,n^))(i-/(x,n' 5 ; 



(15) 



, \T f jStt' .3 n / 1 pi d 9 (T nB -,p+X (y/s) 

\kJ 3 4 ^(V5) p dm> 3 dm> A dm> p 

(i-f(x,W 3 ) ) (i-/(x,n' 4 ) ), (16) 

where the quantities IT (i = 1, .., 5) denote the in-medium momenta of the participating 
baryons. Ilp and Ep stand for the p effective momentum and energy while s = (Hi + 
II2) 2 is the squared invariant energy available in the corresponding baryon-baryon or 
pion-baryon collision. The first sum of the Eq. (15) is that over all baryon-baryon 
collisions while the second sum is that over all pion-baryon collisions occuring during 
the dynamical evolution. An integration over the impact parameter then yields the 
Lorentz invariant differential production cross section 

While in free space the threshold for the elementary production reaction is 4 x 
the nucleon restmass one has to take into account the selfenergies of all participating 
particles for a production process in the medium. The conservation of energy and 
momentum has to be guaranteed, i.e. in baryon-baryon reactions 

p?+p£ = +p£+p£, (is) 

which in terms of effective momenta and effective masses (cf. eqs. (5) and (^j)) reads 

n? + u£{\ pi \,x) + w 2 + ut(\ p 2 \,x) = n£ + u£(\ p 3 \,x) + n£ + 

ua{\ &\, x ) + n% + u»{\ p 5 1, x) + w p + u^{x). (19) 
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With the abbreviation 



A^ = U£(\p 3 \,x) + U£(\p 4 \,x) + U£(\p B \,x) + U^{x) 

-UZQfr \,x)-U£(\p2 U) (20) 

we obtain in shorthand form 

+ n£ = + n£ + n£ + n£ + a**. (21) 

Similar relations are derived for the irN production channels where we assume the pion 
selfenergy to be identically to zero. This assumption might appear questionable |28 



29| , pO] , however, it has been found in a couple of studies that the pion dynamics do 
not indicate strong potentials in the nuclear medium |TT], R3]. 

In order to derive an expression for the differential p-multiplicity we assume, as 
in refs. [|T7|, [H], |3B|, that the differential elementary p production cross section is 
proportional to the phase-space available for the final state in BB reactions: 

d 12 a B B^NNN+p (\/f) f /-v 1 



and 



<5 4 (n^ + - - - rig — — a m ). (22) 



p P P <Pv*b^nn+p (y/s) r 



dm 3 d 3 u 4 d 3 u p — ^ 8jR3(v ^) 

v 



<5 4 (n? + n£ - n£ - n£ - m - A£) (23) 



in case of irB reactions. 

Here, the (^-functions guarantee the energy and momentum conservation and y/s is 
the invariant energy available for the quasi-particles in the initial state. R 4 (\fs) (R 3 ) 
is the 4-body (3-body) phase-space integral [0; it has been included to ensure that 



the differential cross sections are normalized to the total cross section. 



3 Elementary reaction cross sections 



3.1 nN^pX 

Since there are no experimental data on the antiproton production yields in pion in- 
duced reactions close to the reaction threshold [ff5H , we have to address to microscopic 
models in order to obtain closer constraints on this quantity. Here we calculate the cross 
section of the reaction ttN — > pNN within the One-Boson-Exchange model (OBE) us- 
ing the OBE parameters proposed in [^, [37]. Though the OBE results depend on 



the cut-off parameters A* in the formfactors - which leads to an uncertainty in the p 
production cross section by about a factor of 2 according to our investigations - we here 
adopt the results of the Bonn-Julich model 



36 



J for a leading order computation. 
As in Refs. p3L |38[ we will complement the inclusive production of antiprotons at high 
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Figure 1: Processes contributing to antiproton production in nN collisions. 



energies by calculations within the LUND string model |39| (LSM) to include addition- 
ally a multi-hadronic final state with more than two nucleons and an antiproton in the 
final state. 

Fig. |l] illustrates the relevant diagrams for the exclusive reaction ttN — > pNN which 
we will compute in the OBE approach. We take into account both pseudoscalar and 
vector mesons exchanges and express the matrix element as 



T 



nN^pNN 



( . \ \^ ujj/) Tj u{p) Fj(t) 



rrn 



(24) 



where fm-^pN(si) is the amplitude for the process rri — ► pN and the index i stands 
for the exchanged rr, p or cj-meson (cf. Fig. 1), while Fj stands for the corresponding 
formfactor. Here u(p) and u(p') are the spinors of the initial and final nucleons with p, 
p' being their four-momenta, respectively. 

In Eq. (|24D s denotes the square of the colliding energy and si is the square of the 
invariant mass of the pN pair, while t is the four momentum transfer from the initial 
to the final nucleon shown by the lower vertex of the diagram in Fig. [l]. 

The vertices Tj have the following form 



Tjr — g-xNN 75, — gujNN Jfi 



Table 1: The cut-off parameters and coupling constants used in our calculations. 
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A; [GeV] 




fiNN/ 9iNN 


7T 


0.7 


14.4 




P 


1.4 


0.84 


6.1 




1.5 


20 





-p . fpNN . , v 

1 p,p = 9pNN lp + -7T— * °pv K 



2m 



(25) 



where k is the four-momentum of the exchanged p-meson. 

The pNN vertex can be rewritten equivalently using the Gordon's relation for the 
current 



(g P NN + J p nn) lp — 



fpNN 

2m 



(26) 



where m is the nucleon mass, g^NN, 9u>nn, 9 p nn are pseudoscalar and vector coupling 
constants of ttNN, ujNN and pNN interactions, respectively, while J p nn stands for 
the tensor coupling constant. We use the coupling constants from []37|, |4T], [fl] ; f° r 
completeness they are given in Table 1. 

Now the cross section of the reaction ttN — > pNN can be written as 



1 



a 



(s/s-m) 2 



dsx A 1/2 (si, m 2 ,m 2 ) cr^^jv(si) 



lir 2 X(s, m 2 , m 2 ) Jim 2 
where m.j is the mass of the exchanged meson and 

A(X, y, z) = (x - (yV + y/z) 2 ) (x - (y/y - y/zf 

while 

1 



* + A- F 2 (t) 
dt AA| (27) 



t- (t - m 2 ) 



(28) 



t 



± 



2m^ 

2s 



(s + m 2 -m 2 n )(s + m 2 - Sl )T X 1/2 (s,m 2 ,m 2 n )\ 1/2 (s,m 2 ,si) . (29) 



We use the formfactors in the form 



A 2 -m 2 
A 2 -t 



(30) 



with cut-off parameters Aj from |$7], |4lJ (cf. Table 1). 

Note that A„- = 0.7 GeV is taken from j|l] and corresponds to the OBE approach 
of our present paper, while A n = 1.3 GeV is relevant for two-pion exchange corrections 



as considered in Ref. 137 



In Eq. (|27|) the functions At have the following forms: 



— qInn I * 



K = 0Lnn I 2t + 4m 2 |, 



A, 



9 2 p nn (2t + 4m 2 ) + f 2 NN (3t + 5m 2 + — 



+ 8t QpNN f, 



pNN 



(31) 
(32) 
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Figure 2: Total and elastic cross sections for pp interactions as a function of the lab- 
oratory antiproton momentum. The experimental data are from |$5] while the lines 
show the parametrization from the Particle Data Group [(42| . 

In order to calculate cr^i-^pN we address to available experimental data for the 



inverse reaction pp — > nn, np as in Ref. |24j which are related by detailed balance as 

4 si — 4m 2 



a 



TTl— >pp 



Si 



2J t + l X(s u mlml) 
where J, is the spin of the exchanged meson and o^^i is parametrized as [[43|, |4j] 

X^ 2 (s 1 ,m 2 i ,ml) 



(33) 



a 



pp—>1Tl 



si) = Co 



x exp \ -A fsi - K + mi) 2 l 1/2 }> o tot (p lah ) 



(34) 



where C = 9 mb(GeV/c) 1 , A = 4 GeV 1 and a tot {pi ab ) is the total cross section for 
pp interactions shown in Fig. 0. It can be parametrizied by 



a tot (p lab ) =E + Bpf ab + Cln\p lab ) + Dln(p lab ) [mb] 



(35) 
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where pi a b is given in GeV/ c and stands for the antiproton momentum in the laboratory 
system, while the parameters E = 38.4, B = 77.6, n = —0.64, C = 0.26, D = —1.2 are 
taken from the Particle Data Group [^2 . 




Figure 3: The cross sections for the reactions ir+ meson — > pp. The data are extracted 



from the experimental cross sections for the inverse reaction p5| ; the lines show the 
parametrization 



Fig. ^ shows the experimental data for the reaction tt+ meson — > pp calculated via 
detailed balance together with the parametrization (|34|). In the following calculations 
we assume that o-^i^pN is the same for protons and neutrons and does not depend 
upon the charge of the exchanged meson. 

The cross section of the reaction Tt~p — > pnp calculated with Eq. (^7]) is shown in 
Fig. H together with the two available experimental points at higher energy. The solid 
line indicates the total result, while the dashed, dotted and dashed-dotted lines show 
the separate contributions from n, u and p-meson exchanges, respectively. 

The inclusive cross section for the reaction irp — > pX is calculated within the 
framework of the LSM |M| and is shown in Fig. |5| by triangles. Here the dashed line 



indicates the results from the OBE model for the exclusive channel. Fig. | illustrates 
a quite reasonable agreement between the results from two different approaches used 
in our calculations for energies close to threshold. 
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Figure 4: Cross sections for the reactions n~p — > pnp. The experimental data are 
from [35| while the lines indicate our calculations within the One-Boson-Exchange 
model. 

For practical use in transport calculations we parametrize the antiproton production 
cross section in the form 

s \ i /s u 



VnN^pX = a ( 1 — (36) 



with parameters listed in Table. 2. Using the form (p6|) we also parametrize the cross 
section for antiproton production from nucleon-nucleon collisions as calculated within 
the OBE-approach in Ref. ||24j| . Note that the first term of Eq. (|36"D reflects the energy 
dependence of the cross section near the reaction threshold; our final results indicate 
a phase-space dominance for the antiproton production cross sections both from pion 
and nucleon induced reactions which supports the ansatz (23) for the muli-differential 
cross section. It is clearly seen from a comparison of Eq. ( |36| ) with parameters b from 
Table 2. 

Fig. |6| shows our parametrizations of the cross sections for antiproton production 
from np and pp collisions as a function of the excess energy yfs — ^/sq. The solid circles 
indicate the experimental data for the reaction pp — > pX | p5|| . We are not aware of 
inclusive data for p-production by pion induced reactions. Note that the ttN channel 
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Figure 5: The cross sections for the inclusive antiproton production from pion-nucleon 
collisions. The triangles show the results from the string model (LSM). The dashed 
line is our calculation for the exclusive reaction within the OBE approach while the 
solid line indicates the parametrization (|36|). 

Table 2: The parameters in the approximation fl36|). 



Reaction 


so 


a [mb] 


b 


c 


irN -^pX 


9 m 2 


1 


2.31 


2.3 


NN — > pX 


16 m 2 


0.12 


3.5 


2.7 



above threshold is much larger than the pp channel since it increases with 3-body phase 
space compared to 4-body phase space for pp. This indicates that the contribution 
from secondary pion induced reactions to antiproton production in proton-nucleus and 
heavy-ion collisions might become important or even of leading order. 

In the following calculations we assume that the cross sections are the same for the 
reactions with protons, neutrons and A- resonances and for ir + , n~ and 7r°- mesons. 
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Figure 6: Cross sections of inclusive antiproton production from up and pp collisions. 
The circles show the experimental data for the reaction pp — > pX |35] while the lines 
are the parametrizations (|36|) of our calculations. 



4 Comparison to experimental data 



Following the earlier studies on p production |22], we now attempt to extract in- 



formation on the antiproton potentials in the medium by comparing our transport 
calculations with the available data, i.e. we try to fix the coupling constants gf , g% 
empirically first starting with g% — g% — 0. 

4.1 Proton- nucleus collisions 

Detailed experimental studies on antiproton production in p + A collisons at beam 
energies below the free NN — > NNNp reaction threshold were performed at KEK ]TT| 
45| . Antiprotons with momenta of 1. - 2.5 GeV/c were detected at an emission angle of 



5.1° in the laboratory and at beam energies of 3.5, 4, 5 and 12 GeV. The experiments 
show a very high production cross section at the incident kinetic energy of 3.5 GeV, 
which is substantially below the free iViV threshold of 5.6 GeV. 

Here we perform the analysis of antiproton production within the transport model 
described in Section 2. For the propagation of antiprotons in nuclear matter we account 
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for both the absorption and elastic scattering. We assume that the elastic cross section 
is 37% of the total cross section; this simplifying assumption is shown in Fig. ^| by the 
solid line which reasonably fits the experimental data. 

The antiprotons - once produced - strongly interact with the surrounding nuclear 
matter due to the large pN cross section and a large fraction annihilates. Fig. [7| 
illustrates the nuclear transparency to antiprotons calculated for different targets and 
defined as the ratio of antiprotons detected asymtotically for t — > oo to the total 
number of antiprotons actually produced in baryon-baryon and pion-baryon collisions. 

The histograms in Fig. ^ show our results from the transport model while the solid 



lines indicate the calculations within the Glauber formalism described in |3q| . The 
Glauber model accounts for the interaction of the initial protons as well as for the 
annihilation of the antiprotons. However, the rescattering in the nuclear medium is 
not included in the present version of the Glauber model. Thus the large difference at 
antiproton momenta p > 2 GeV/c between the Glauber calculations and the results 
from the BUU approach illustrated in Fig. |7| is due to rescattering effects, which are 
most important for heavy nuclei. We find that both calculations agree reasonably well 
for 12 C and 64 Cit for which experimental data are available. Both models predict that 
for heavy nuclei in p + A reactions about 90% of the produced antiprotons annihilate 
in the target nucleus. 

Following the investigations from |22|, |23| we study the influence of the antiproton 
self-energy on the differential p-production. In the numerical calculations we adopt the 



limit g% = as in |Z3|| and vary only the scalar coupling g& in Eq. (4) in order to obtain 
the best fit to the experimental data. 

Fig. |8] shows the invariant differential cross section for antiproton production in 
p + 12 C collisions at a beam energy of 3.5 GeV and an emission angle of 5° in the lab- 
oratory. The experimental data are taken from |^5| while the histograms indicate our 
results calculated with different antiproton potentials (in MeV at saturation density 
Po = 0.16 fm~ 3 ). It is clearly seen that without in-medium effects (U = 0) we signifi- 
cantly underestimate the production cross section. The agreement becomes better for 
an attractive antiproton potential of —125 ±25 MeV at normal nuclear matter density. 
Figs. |5|JT(]|, furthermore, show the experimental data together with our calculations 
for the reaction p + C — > pX at bombarding energies of 4 and 5 GeV, respectively. 
Note that at higher energy the data are no longer that sensitive to a variation of the 
antiproton potential. 

The antiproton spectra from p + Cu collisions at bombarding energies of 3.5, 4 and 
5 GeV are shown in Fig. 11. The dashed histograms are the results for U — MeV, 
while the solid histograms indicate our calculations with an antiproton potential of 
U = —100 MeV at po which almost perfectly fits the experimental data at all energies. 

As anticipated before, the nN production channel dominates in p + A reactions at 
least up to bombarding energies of 6 GeV (cf. Table 3). This finding is fully in line 
with our studies on K + , K~ , p, u and production for proton-nucleus reactions p8| 
and results from the fact a) that the Fermi motion can be exploited twice in secondary 
reactions and b) that the secondary tcN collisions have much larger production cross 
section close to threshold (cf. Fig. 6). 

In summary, we get a reasonable fit to the experimental data on antiproton produc- 
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Figure 7: Nuclear transparency for antiprotons produced in p + A reactions. The 
histograms show calculations within the transport model while the solid lines indicate 
the results from the Glauber approach in Ref. [BH . 



Table 3: Cross section a [nb] for antiproton production in p + 12 C collisions at bom- 
barding energy T p calculated for different p potentials U [MeV]. 



T p [GeV] 


Reaction 
channel 


U [MeV] 





-50 


-100 


-150 


-200 


5.0 


nN — > p 
NN — > p 


362 
8.4 


563 
20.5 


824 
41.7 


1210 
78.5 


1670 
138 


6.0 


ttN — > p 
NN -+p 


2810 
250 


3784 
384 


4895 
569 


6361 
816 


8076 
1143 



tion from proton-nucleus collisions at subthreshold energies with an attractive antipro- 
ton potential U ~ — 125 ± 25 MeV at normal nuclear matter density and at antiproton 
momenta from 1 to 2.5 GeV/c with respect to the nuclear matter at rest. 
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p + 12 C^p + X at 3.5 GeV 
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Figure 8: Antiproton spectrum from p + C collisions at a beam energy of 3.5 GeV. The 
experimental data are from while the histograms indicate our calculations with 
different antiproton potentials U in MeV. 

4.2 Nucleus- nucleus collisions 

Experimental studies on antiproton production from heavy-ion collisions at energies 
below the free NN threshold have been performed at the BEVALAC JBj, [|7| and 



at GSI Darmstadt |L2 



Whereas in proton-nucleus collisions the antiprotons have quite large momenta rel- 
ative to the nuclear matter, in heavy-ion reactions the antiprotons have small momenta 
in the nucleus-nucleus center of mass and are comoving with the expanding 'fireball'. 
Moreover, heavy-ion reactions probe much higher baryon densities, which influences 
both the antiproton selfenergy and their final state interactions, i.e. reabsorption and 
elastic rescattering. 



Fig. [12] shows the antiproton spectra from Si + Si collisions at a bombarding energy 
of 2 A GeV. The experimental data are from |1(], [|7| while the histograms indicate our 
calculations for different values of the antiproton potential U (in MeV at po). In 
Table 4 we illustrate the relative contribution to antiproton production from NN, AN 
and ttN reaction channels for a potential U = —150 MeV at po- 111 order to investigate 
the variation of the cross section with the impact parameter b we also show in Table 4 
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p + 12 C^p + X at 4.0 GeV 
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Figure 9: Antiproton spectra from p + C collisions at a beam energy of 4 GeV at 
Oiab = 0°. The experimental data are from |45j while the histograms indicate our 
calculations with different antiproton potentials U in MeV at density po- 



the differential antiproton multiplicity multiplied by the factor 2nb. It becomes clear 
that the dominant contribution to p-production also for A + A collisions stems from the 
secondary pion induced reactions as in case of p + A collisions. Moreover, the dynamics 
of antiproton production and propagation show no sizeable difference in the channel 
decomposition as a function of the centrality of the collision. 



Fig. [13] shows the antiproton spectra from Ni + Ni collisions at 1.85 A GeV. The 
experimental data are taken from [TB|, and can be reasonably reproduced by our 
calculations with U ~ —125 ± 25 MeV at po whereas the data are underestimated by 
about 2 orders of magnitude if no p potential is included. Similar statements also hold 
for the antiproton spectra from Ne + NaF collisions at 1.94 A GeV as shown in Fig. |T4 
in comparison to the preliminary experimental data from . 



4.3 Dispersion approach for the antiproton potential 

The antiproton potential extracted above empirically from experimental data on p- 
production from p + A and A + A reactions now can be compared with the previous 
estimates from p2[ Q as well as with the optical potential from a dispersion relation 
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p + 12 C^p + X at 5.0 GeV 
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Figure 10: Antiproton spectra from p + C collisions at a beam energy of 5 GeV. The 
experimental data are from while the histograms indicate our calculations with 
different antiproton potentials U in MeV at po- 

Table 4: The differential antiproton multiplicity for Si + Si collisions at 2 A GeV 
multiplied by the factor 2-7T& for an antiproton potential of U = —150 MeV at po- 
The decomposition is performed for different impact parameter b and nucleon-nucleon, 
A-nucleon and pion-nucleon reaction channels. 



b [fm] 


da/db [nb/fm] 


NN 


AiV 




1 


0.64 


7.09 


22.0 


2 


1.1 


9.5 


28.9 


3 


1.0 


7.3 


24.5 


4 


0.67 


3.3 


4.7 


5 


0.027 


0.63 


1.7 



approach. In the latter approach the real part of the Schrodinger equivalent potential 
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Figure 11: Antiproton spectra from p + Cu collisions at beam energies of 3.5, 4 and 
5 GeV. The experimental data are from [45] while the solid histograms show our cal- 



culations with an antiproton potential U = —100 MeV at po; the dashed histograms 
indicate our results for U = 0. 

in the low density approximation is calculated by means of the dispersion relation as 

U(E, PB ) = -E* P r dE ' f^g^w (37) 
V H ' 2tt Jm E' (E-E') K J K ' 

where P stands for the principal value of the integral, a ann is the annihilation cross 
section while ps is the baryon density. 

The solid line in Fig. |T5] shows the antiproton potential calculated by ( B7| ) using the 
parametrization of the annihilation cross section from [|L4}| . We note, that the integral 
(|3~7D depends sensitively on the parametrization of the annihilation cross section at 



high energy and that the present result using Ref. [Q should be regarded cautiously. 
Fig. [15], furthermore, also shows our results extracted from p + A and A + A reactions 
by the rectangles, which are in a rough agreement with (^) at low momentum. 
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Figure 12: Antiproton spectra from Si + Si collisions at 2.0 A GeV. The experimental 
data are from [f|7]] while the histograms show our calculations with different antiproton 
potentials (at po). 



5 Summary 



In this study we have evaluated the differential cross section for p production for 
proton-nucleus and nucleus-nucleus reactions in the subthreshold regime by consider- 
ing on-shell baryon-baryon production channels involving nucleons and A's with their 
in-medium quasi-particle properties and pion-baryon collisions while treating the p 
propagation and annihilation explicitly The quasi-particle properties of the nucleons 
are fixed within our approach by the nuclear saturation properties, the proton-nucleus 
empirical potential as well as Dirac-Brueckner calculations at higher density as in 
|f22| [2~3| , |32|| . Whereas previous studies have neglected the tiN production channel and 
used extrapolations to threshold of the elementary process p+p — > p + X, our approach 
is based on the elementary productions cross section calculated within the OBE ap- 
proach and within the Lund-string-model at higher energy. These cross sections are in 
agreement with 3- and 4-body phase-space contrary to the extrapolation from Batko 
et al. |19| - used so far in almost all studies - which overestimates the pp channel 
substantially close to threshold. 

We have performed a systematic study of p-nucleus and nucleus-nucleus collisions 
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Figure 13: Antiproton spectra from Ni + Ni collisions at 1.85 A GeV at Q\ a b = 0°. The 
experimental data are from fll3| , pf] while the histograms show our calculations with 
different antiproton potentials U at density po — 0.16 fm~ 3 . 



in a broad kinematical regime and compared our numerical results to the data from 
KEK [|lTJ and GSI []12" , 13|. We find a consistent description of all data employing an 
attractive potential for the antiprotons of about -100 to -150 MeV at po which is roughly 
in line with a dispersive potential extracted from the dominant imaginary part of the 
antiproton selfenergy due to annihilation. Since the dispersive model itself is known to 
be valid only at small baryon density and the medium corrections to the annihilation 
cross section are practically unknown, the agreement between our empirical analysis 
and the dispersion approach is not considered to be a strong argument. Whereas our 

@the 
•125±25 



novel analysis roughly yields the same p potential at low p momenta as in Ref. 
antiproton potential should be more attractive for momenta above 1 GeV/c 
MeV) than assumed before. 

Our analysis at subthreshold energies, on the other hand, suggests lower antiproton 
potentials as anticipated from studies at AGS energies |50, 51] where the total antipro- 



ton production cross section is no longer that sensitive to its value near the reaction 
threshold. Here a p potential of ~-250 MeV at po is proposed in Refs. p0| , |5T| for an 
antiproton at rest with respect to the nuclear matter, whereas an antiproton potential 
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Figure 14: Antiproton spectra from Ne + NaF collisions at 1.94 A GeV. The circles 
show the preliminary results from ||9| while the histograms indicate our calculations 
for different antiproton potentials (at po). 



of ~-160 MeV is reported for a p momentum of 1 GeV/c |5TJ . 

Futhermore, independent from the above uncertainties it has become clear that 
the p production at subthreshold energies is dominated by secondary pion induced 
reactions as in case of kaons and antikaons f25| , |2~6fl since for the new elementary cross 
sections the baryon-baryon channel is suppressed substantially 



The authors acknowledge valuable discussions with A. Gillitzer and J. Chiba through- 
out this study and for providing us with their reanalyzed experimental data prior to 
publication. 
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